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FOREWORD 


This  tr.hnii-.il  memorandum  is  .1  copy  i>t  an  article  submitted  for 
puhlir.it  ion  to  the  Institute  for  Electrical  and  Electronics  Engineers 
(IEEE)  Transactions  on  Antennas  and  Propagation.  Some  of  the  results 
presented  herein  plus  additional  results  on  geometry  deformations  due 
to  atmnspher tc  refraction  has  been  submitted  to  the  Union  Rad io-Se lent i f lque 
Internationale  (t'RSl)  triennial  Symposium  on  Electromagnetic  Wave  Theory, 
which  is  to  be  held  in  conjunction  with  the  1FF.E  Internationale  Symposium 
on  Antennas  and  Propagat ion  and  the  annual  technical  meeting  of  the  U.S. 
National  Committee  of  URSI.  The  Electromagnetic  Wave  Theory  Sumposlum 
will  h«-  held  on  2 0-24  .June  1977  at  Stanford  University,  Palo  Alto, 
California.  The  work  presented  in  this  paper  has  been  supported  by  a 
contract  with  the  Naval  Air  Systems  Command  (AIR-170C). 

This  document  has  been  prepared  as  an  interim  presentation  of  timely 
information.  Although  care  has  been  taken  in  the  preparation  of  the 
technical  material  presented,  tin-  results  herein  are  to  be  considered  as 
prelimlnarv  in  nature. 


PAH.  B.  HOMER 

Head,  Weapons  Systems  Analysis  Division 
Systems  Development  Department 
S  .January  1976 
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ABSTRACT 


The  derivation  of  a  transform.it  Ion  to  compensate  for  refract lvlty 
In  the  earth's  atmosphere  Is  obtained  through  the  use  of  ray  tracing 
and  differential  geometry.  The  effective  earth's  radius,  which  is  the 
classical  method  to  compensate  for  the  ray  curvature  as  an  electro¬ 
magnetic  wave  propagates  through  the  earth's  atmosphere,  is  also  derived 
and  shown  to  be  invalid  since  it  does  not  serve  its  Intended  purpose  — 
namely,  to  conformally  transform  curved  rays  into  straight  lines.  A 
brief  rev i ew  of  a t mospher i c  radio  refractivity  and  rav  tracing  is  pre¬ 
sented  and  the  differences  between  the  refraction  compensation  transforma¬ 
tion  and  the  effective  earth's  radius  are  demonstrated. 
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altitudes  of  the  ray's  end  points 

refractlvlty  constants 

absorption  Index  of  refraction 

effective  spacing  between  neighboring  rays 

complex  Index  of  refraction 

refract Ivlty 

refractlvlty  at  the  earth’s  surface 
phase  Index  of  refraction 
total  atmospheric  pressure 
carbon  dioxide  partial  pressure 
atmospheric  pressure  of  dry  air 
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equatorial  radius  of  the  earth 
polar  radius  of  the  earth 
radius  of  the  earth 
effective  radius 
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rudiuH  of  curvature  of  surface  of  constant  refract Ivlty 
transformed  radius 

distances  from  earth's  center  to  the  ray's  end  points 
atmospheric  temperature 
ray  length 

refraction  compensation  transformation  constant 
gradient  operator 
dielectric  constant 
ray  factor 

magnetic  permeability 

angle  ray  makes  with  surface  of  constant  refractivity 

.yigles  that  the  rav  makes  with  the  surface  of  constant 
refractivity  at  the  ray's  end  points 

latitude  on  the  earth’s  surface 

radius  of  curvature  of  the  ray 

earth-center  angle 

transformed  earth-center  angle 
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INTKODUC1 ION 


From  the  Fermat's  principle  of  least  action,  electromagnetic 
radiation  (l.e.,  photons)  always  propagate  along  a  geodesic  in  free 
space  or  through  any  homogeneous  medium.  However,  if  the  radiation  is 
propagating  through  an  inhomogeneous  medium  (e.g.,  the  earth's  atmos¬ 
phere  ) ,  gradients  in  the  speed  of  light  in  the  medium  cause  the  path  of 
the  photons,  or  ray,  to  deviate  f rom  a  geodesic,  or  straight  line.  If 
the  speed  of  light  In  the  medium  Is  discontinuous,  the  ray  path  will  not 
he  differentiable  (e.g.,  at  reflection  points);  however,  for  the  earth's 
atmosphere  it  Is  reasonable  to  assume  that  the  speed  of  light  is  a 
differentiable  function  of  location.  Then,  the  rays  are  twice  differ¬ 
entiable  curves  through  the  atmosphere  and  hence  have  a  curvature 
defined  along  the  ray  path. 

This  lurvaturc  of  the  ray  paths  due  to  inhomogenelt les  in  the 
atmosphere  causes  considerable  difficulties  In  attempting  to  solve 
various  geometric  problems.  Although  rav  tracing  is  one  solution,  it  is 
usually  much  too  complex  for  most  applications.  Another  technique  has 
been  the  development  of  a  t ransform.it  Ion  that  maps  the  curved  rays  into 
straight  lines  and  preserves  with  sufficient  accuracy  various  geometric 
parameters  such  that  the  transformed  geometry  can  he  used  to  solve 
certain  problems.  The  most  cotmnon  method  used  to  compensate  for  the  rav 
curvature  In  the  earth's  atmosphere  is  to  use  an  effective  earth's 
radius,  which  Is  usually  larger  than  the  true  earth's  radius,  that 
supposedly  conformally  transforms  the  curved  rays  into  straight  lines. 
The  treatment  of  rav  curvature  and  the  effective  earth's  radius  in  the 
standards  of  the  literature  has  not  been  sufficiently  rigorous  due  to 
either  unnecessary  or  tacit  assumptions  |1)  -  [4].  This  lncomp leteness 
leads  to  incorrect  results  due  to  improper  use  of  the  effective  earth's 
radius.  In  addition,  the  effective  earth's  radius  Is  not  a  conformal 
t ransformat Ion,  and  It  does  not  transform  the  curved  ravs  Into  straight 
lines.  Instead,  the  curved  rays  are  transformed  Into  ravs  with  a 
greater  radius  of  curvature. 

In  this  paper,  the  basic  theory  of  atmospheric  radio  refractlvity 
Is  presented.  Although  the  concepts  of  rav  curvature  and  effective 
earth's  radius  are  applicable  for  all  frequencies,  this  paper  con¬ 
centrates  on  the  propagation  of  microwave  radiation.  After  a  brief 
discussion  of  geometric  optics,  the  effective  earth's  radius  is  derived 
and  son.  of  the  associated  errors  are  discussed.  Finally,  a  trans¬ 
formation  that  compensates  for  the  ray  curvature  Is  derived  and 
compared  with  the  effective  earth's  radius. 
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REFRACT  I  Vm 


The  index  of  retraction  of  a  medium  is  defined  to  be 

m  •  CU  ■  n-lk,  (1) 

where  •  is  the  dielectric  constant  and  u  is  the  magnetic  permeability  of 
t tie  medium.  The  real  part  of  the  index  of  refraction,  n,  called  the 
phase  index,  is  the  ratio  of  the  speed  of  light  in  a  vacuum  to  that  In 
tiie  medium,  while  the  imaginary  part  of  the  index  of  refraction  Is 
the  absorption  index,  k.  For  the  remainder  of  this  report,  the  absorp¬ 
tion  index  will  be  Ignored.  For  electromagnet ic  radiation  below  30  GHz 
in  the  earth’s  atmosphere,  the  phase  index  does  not  exceed  unity  by  more 

than  4  SO  x  10  h.  so  a  parameter  called  refract ivlty  is  defined  by 

N  -  (n-1)  x  106.  (2) 

For  frequencies  below  10  GHz,  the  refractivlty  is  given  by  Oebye 
theory  for  polar  and  non-polar  molecules  as: 


:1  T  *  K2  T  +  K1  ?2  *  K4  T 


(3) 


where  P  is  the  atmospheric  pressure  of  drv  air,  P  is  the  partial 

d  c 

pressure  of  carbon  dioxide,  e  Is  the  partial  pressure  of  water  vapor, 
and  T  is  the  atmospheric  temperature  Ill.  However,  since  In  the  tropo¬ 
sphere  (t.o.,  for  altitudes  less  than  10-11  kilometers)  the  concentration 
of  carbon  dioxide  bv  molecular  weight  is  approximately  0.03  percent  and 

K./K,  1/3,  it  is  customary  to  combine  the  pressure  of  dry  air,  P., 

s  l  a 

with  the  partial  pressure  of  carbon  dioxide,  P^ ,  resulting  In 


N  -  Ki  r  +  k2  ? +  *3  ^  • 


(4) 


where  P(  is  the  total  atmospheric  pressure.  F.xper imental  data  results 
in  the  following  values  for  the  constants  | 1 1 : 

K/mb 


K  J  •  77.607  +0.13 
K,  •  71.6  +8.5 

-  (3.747  +  0.031)  x  lo' 


K/mb 

K:/mb 
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Another  common  expression  for  the  ref  ract  lvity ,  which  is  sufficiently 
accurate  for  many  applications,  is  (1) 


N 


e 


2 

T 


where 

Ks  -  3.71  x  10s  IC/mb 


(5) 


The  constants  in  (5)  are  considered  to  be  accurate  to  within  0.5  percent 
in  N  for  frequencies  up  to  10  (’.Hi  and  normally  encountered  ranges  of 
pressure,  temperature  and  humidity  II).  The  first  term  of  (5)  is  usually 
referred  to  as  the  "dry"  term  and  the  second  term  is  referred  Co  as  the 
"wet"  term.  As  can  he  seen  from  (*)  and  (5),  the  propagation  of  micro- 
wave  energy  through  the  earth’s  atmosphere  is  affected  primarily  by 
gradients  in  the  wet  term. 


ATMOSPHERIC  REFRACT  I VITY  MODELS 


Although  t he  expressions  for  refractivity  as  a  function  of  pressure, 
temperature,  and  water  vapor  partial  pressure  are  essential  if  one 
desires  to  evaluate  electromagnetic  wave  propagation  under  particular 
atmosphor Ic  conditions,  it  is  frequently  convenient  to  use  a  model  that 
represents  the  mean  refractivity  profile  for  a  given  location  or  region. 
Several  models  in  Common  use  today  are  discussed  below.  It  should  he 
mentioned,  though,  that  the  deviation  from  the  mean  may  he  quite  signifi¬ 
cant  and  the  profile  that  is  referred  to  as  "normal"  may  occur  only  a 
small  percentage  of  the  tins'. 

LINEAR  ATMOSPHERE 

The  linear  model  assumes  that  the  mean  refractivity  varies  linearly 
with  the  altitude;  t.e., 

N  *  N  ♦  h  4r  ,  (€>) 

o  dh 

where  N  is  the  mean  refractivity  at  zero  altitude,  and  h  is  the 

o 

altitude.  This  model  can  only  he  considered  valid  for  low  altitudes 
(i.e.,  less  than  1  kilometer).  Standard  conditions  are  considered  to  he 

N  -  301.0  and  39.25  km'1, 

o  dh 
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EXPONENT !  A!  ATMOS  PH  F.k  K 

In  this  model,  the  me. in  ref  met  ivlt  v  varies  exponentially  with  the 
n  1 1 1 1  udo ;  l .  e . , 

-c  h 

N  -  N  e  °  ,  (7) 

o 


where  c  is  a  constant.  This  model  gives  f.iirlv  close  agreement  with 

o 

the  mean  refract ivlty  for  altitudes  less  than  3  kilometers  and  is  still 
moderately  close  for  much  higher  altitudes.  The  exponential  atmosphere 
has  been  adopted  bv  the  National  Bureau  of  Standards  with  values  of 

and  e  given  in  Table  1  |1).  With  these  values  the  exponential  model  Is 

o 

referred  to  as  the  Central  Radio  Propagation  Laboratory  (CRPL)  Exponential 
Reference  Atmosphere  (1958). 

OTHER  MEAN  ATMOSPHERES 

Other  model  atmospheres  have  been  developed  to  represent  the  mean 
refractivlty  as  a  function  of  altitude.  One  such  model  uses  various 
combinations  of  linear  and  exponential  atmospheres  for  different 
altitude  ranges.  Another  model  uses  a  b 1  exponent i a  1  atmosphere,  in 
which  the  wet  and  dry  terms  of  ( *> )  are  assumed  to  varv  exponentially  but 
with  different  scale  height  constants.  These  model  atmospheres  can  be 
found  in  the  Reference  [1|. 

NON-STANDARD  ATMOSPHERES 

The  use  of  mean  and  standard  atm.  spheres  has  created  much  erroneous 
data  and  misconceptions  concerning  the  propagation  of  microwave  radia¬ 
tion  in  the  earth's  atmosphere  since  many  operational  personnel  have 
treated  the  standard  atmosphere  as  the  most  likelv  one.  Non-standard 
atmospheres  which  involve  various  anomalous  effects  are  quite  common  in 
various  parts  of  the  world  and  have  a  significant  effect  upon  the 
propagation  of  electromagnetic  waves  ( S ) .  A  studv  of  these  non-standard 
atmospheres  is  bevond  the  scope  of  this  paper;  however,  certain  micro- 
wave  propagation  conditions  can  be  classified: 


o.n 
0.5 
1.0, 

where  P  is  the  rav  factor  defined  by 
rn  dn 

‘  "  n  dr”  ’ 
n 


Sub  refraction:  0.0  <  T 

Normal  refraction:  -0.5  <  T  < 

Superrefraction:  -1.0  <  T  <  - 

Trapping  (ducting):  P  <  - 


(8) 


9 


NVK:  I'M  120H 

lAHi.K  t.  Table  of  Constants  for  the  (31  PL  Exponential 
Radio  Refractlvlty  Atmosphere  (1958)  (1). 
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where  n  In  the  phase  Index  of  refraction,  which  Is  assumed  to  be  a 

differentiable  1  unction  of  r  ,  and  r  Is  the  radius  of  curvature  of  the 

n  n 

surface  of  constant  ref ract ivlty.  For  a  spherically  stratified  atmos¬ 
phere  r  »r,  tin-  distance  from  the  earth's  center.  It  should  he  noted 
n 

that  this  classification  distinguishes  between  superrefraction  and 
trapping.  This  distinction  is  not  ;>  1  wavs  made  in  the  literature; 
however,  due  to  tin*  significant  difference  in  propagation  characteristics, 
it  is  felt  that  the  distinction  should  he  made. 


('.Ft  'METRIC  OPTICS  (RAY  TRAC  INC.) 


Since  the  concepts  of  geometric  optics  are  used  in  developing  the 
theory  of  the  effective  earth’s  radius,  it  is  useful  to  evaluate  when 
geometric  optics,  or  rav  tracing,  is  valid  and  to  give  some  expressions 
for  ray  tracing  in  the  earth's  atmosphere.  For  a  more  detailed  treat¬ 
ment  of  the  subject,  the  reader  is  referred  to  Stavroudis  [fe],  an  excellent 
mnthem.it  lea  1  treatment  of  rav  theory.  Rav  tracing  is  a  valid  approxima¬ 
tion  to  the  propagation  of  electromagnetic  waves  through  a  medium  if  the 
phase  index  of  refraction  does  not  change  s  1  gn leant  1  v  within  a  distance 

2 

comparable  to  the  wavelength  of  the  signal  (i.e.,  |Vp  <<  2~n  X)  and  if 
the  relative  change  in  the  effective  spacing,  I.  between  neighboring 
rays  does  not  change  significantly  with  respect  to  the  wavelength  of  the 
signal  (l.e.,  |  VI.  I  •-<  2ttnl./X).  Hoth  of  these  assumptions  are  valid  for 
non-trapping  atmospheric  conditions.  11  trapping  occurs,  special  care 
must  he  iken  to  avoid  caustics  and  various  diffrartion  phenomenon. 

The  assumption  is  made  for  the  remainder  of  this  paper  that  the 
phase  index  of  refraction  is  a  differentiable  function  that  depends  upon 
altitude  only.  The  assumption  that  n  depends  upon  r  onlv  (i.e.,  the 
atmosphere  is  spherically  stratified,  or  horizontally  homogeneous)  is  a 
fatrlv  conventional  one;  although  this  assumption  is  not  necessarv  for 
the  development  of  r  iv  tracing  or  the  effective  earth's  radius,  it 
simplifies  the  mathematics  somewhat.  The  assumption  that  n  is  a 
differentiable  function  is  not  too  stringent  a  condition;  it  mav  slmplv 
require  some  smoothing  of  experimental  data. 

Us  fng  Snel 1 ' s  l  aw, 

nr  cos  $  »  constant,  (9) 

n 

and  some  simple  differential  geometry,  the  above  assumptions  quickly 
give  rise  to  the  following  relations  for  a  rav  propagating  through  a 
horizontally  homogeneous  atmosphere: 


11 
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dr  r  ran  $ 

di  — T*r~ 


I  ! 
d  I  '  •  ‘ 


(11) 


dr 

dO 


r  ran  t , 


<  1  2 ) 


dx  _  _ n _  ,  (13) 

dr  s  1  n  $ 


„  i  -r  cos  <> 

c  ■  -  - - 

P  r 


(14) 


where  r  i  •.  the  distance  from  the  earth's  center,  t  Is  the  angle  between 
the  rnv  and  the  surface  of  constant  ref ract 1 vlty ,  0  is  the  earth-cent er 
angle,  x  is  the  distance  along  the  ray  between  the  ray's  end  points,  C 
is  the  ray  curvature,  ,  is  the  radius  of  curvature  of  the  ray  and  T  is 
the  ray  factor. 

It  should  be  noted  that  for  this  paper,  the  actual  value  of  the 
earth's  radius  1  not  of  great  importance.  This  discussion  is  concerned 
nnlv  with  ravs  propagating  between  two  points  in  the  earth's  atmosphere 
and  the  interaction  of  the  electromagnet lc  wave  with  the  earth's  surface 
is  not  considered.  Therefore,  all  of  the  analysis  herein  is  valid  for 
an  atmosphpre  whose  surfaces  of  constant  refractivlty  are  concentric 
spheres  and  whose  radii  may  have  nothing  whatsoever  to  do  with  the 
radius  of  the  earth.  However,  for  reference  and  for  numerical  data 
presented  near  the  end  of  this  paper,  the  appendix  contains  a  brief 
discussion  of  the  radius  of  the  earth. 


effective  earth's  radius 


The  classical  approach  to  compensate  for  the  curvature  of  rav  paths 
duo  to  atmospheric  refraction  is  to  assume  that  the  earth's  radius  is 
somewhat  larger  than  it  actually  is.  This  method  was  advanced  by  Schelleng, 
Burrows  and  Ferrell  (7)  and  has  been  used  by  many  authors  since.  The 
concept  is  to  absorb  the  curvature  of  the  ray  into  the  curvature  of  the 
effective  earth  so  that  the  relative  curvature  remains  the  same.  If  r 
is  the  distance  from  a  point  on  the  ray  to  the  center  of  the  earth  and  p 
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Is  the  radius  of  curvature  of  the  ray  at  that  point,  the  effective 

radius  r,  is  defined  to  he 
k 

-  (r  1  -  p  1 )  1  .  (H) 


Since  anv  twice  differentiable  curve  has  a  radius  of  curvature  associated 
with  each  point  on  the  curve,  this  technique  applies  to  any  realistic 
rav  path  (i.e.,  if  n  is  differentiable).  The  negative  sign  results  from 
the  convention  that  a  downward  curving  rav  has  positive  curvature  and 
an  upward  curving  rav  has  negative  curvature.  This  results  in  an 
effective  radius  of 


j 

k  i+r  cos  i 


(16) 


Fquatlon  (lh)  Is  applicable  for  any  atmospher lc  condition  in  which 
rav  tracing  is  valid.  However,  both  T  and  $  are  functions  of  the 
position  of  the  point  along  the  rav  path;  therefore,  V  cos  4  varies  along 
the  ray  path.  For  this  reason,  the  effective  radius  is  useful  if  and 
onlv  If  the  variation  in  I'  cos  $  along  the  ray  path  is  small  and  can  be 
ignored.  However,  it  should  be  pointed  out  that  at  low  altitudes  for 
exponential  and  linear  atmospheres ,  r^  decreases  with  increasing 
alt  it ude. 

Manv  authors  have  treated  the  effective  earth's  radius  using 
assumptions  that  are  not  always  explicitly  stated.  The  effective 

earth's  radius,  r,'  Is  defined  as  the  effective  radius  evaluated  at 

the  earth's  surface  with  *  •  0,  so 


where  r  is  the  radius  of  the  earth.  It  is  then  claimed  that  the 

e 

effective  radius  that  should  be  used  Is 


(18) 


since  r  *  r  ♦  h.  It  should  be  noted  that  there  is  significant  dif- 
•  (h) 

ference  between  r^  and  r^  .  For  normal  propagation,  T  ~  -1/4,  so 

(0)  4 

rfc  -  ?  r^  ,  the  so-called  4/3-earth  model.  If  normal  propagation  is 
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assumed,  tor  small  look  angles  ($  0)  differences  of  ))X  exist  between 

the  altitude  difference  of  the  ray's  end  points  calculated  using  the 
effective  earth's  radius  and  that  using  the  effective  radius;  for  a  very 
high  look  angle  ($  v/2),  the  difference  Is  approximately  zero.  Also, 

due  to  the  high  chances  of  aubref ract Ion,  superref ract Ion,  and  trapping 
propagation  conditions  In  certain  areas  of  the  world,  other  values  of  I 
mav  he  necessary.  Depending  upon  the  atmospher lc  conditions,  the  actual 
differences  might  be  less  or  greater  than  those  above. 

It  is  often  assumed  that  the  effective  earth's  model  is  applicable 
only  when  the  linear  atmosphere  Is  valid.  However  the  only  requirement 
for  appl lcablllty  Is  that  T  be  relatively  uniform  over  the  ray  path 
between  the  transmitter  and  the  receiver,  and  the  model  of  at mospher 1 c 
refractlvity  Is  Immaterial. 

I'slng  the  effective  earth's  radius.  It  Is  assumed  that  all  trans¬ 
formed  rays  can  be  represented  as  straight  lines  over  the  surface  of 
the  effective  earth  and  that  the  t ransf ormat ion  Is  conformal  (l.e., 
angles  are  preserved).  However,  If  one  examines  the  transformed  ravs. 

It  becomes  readily  apparent  that  thev  are  not  straight  lines,  and  the 
t ransforraat Ion  is  not  conformal.  Kor  example,  since  the  effective 
earth's  radius  transformation  does  not  alter  the  earth's  center  angle 
between  the  transmitter  and  the  receiver,  it  is  Impossible  to  draw  a 
triangle  that  also  preserves  the  angles,  and  between  the  rav  and 

the  surfaces  of  constant  refract lvltv  at  the  transmitter  and  receiver, 
respectively;  therefore,  if  the  transformed  ravs  were  straight,  the 
t ransforraat Ion  would  not  be  conformal.  Also,  a  detailed  study  of  the 
transformed  rays  show  that  although  their  curvature  has  decreased,  it  is 
not  zero  (l.e.,  it  is  not  a  straight  line).  In  the  following  section,  a 
t rans format  Ion  that  maps  the  curved  ravs  into  straight  lines  Is  derived 
and  the  differences  between  that  transformation  and  both  the  effective 
earth's  radius  anil  the  effective  radius  t ransf ormat Ions  are  demonstrated. 


KKKRACTIt'N  OOMVKNSAT I  ON  TRANSFORMATION 


The  curvature  of  a  twice  differentiable  curve  In  spherical  coordinate 
Is  defined  to  he 


C 


j  i 

r“  ♦  2r‘  -  r  r 


.  2  .2,3/2 

(r  ♦  r  ) 


(19) 


where  f  Implies  differentiation  with  respect  to  d.  Requiring  that  the 
angle  >  be  preserved  all  along  the  ray,  and  using  differential  geometry, 
the  condition  that  the  transformed  rav  is  a  straight  line  (l.e.,  the 
curvature  is  zero)  Is  satisfied  by 
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r 


♦  2r* 


(20) 


which  Is  satisfied  If  and  only  If 


dr  r 

~  m  ~  <I«T>  (21a) 

dr  r 


d8  - 

<-r~)  *  l  .  (21b) 


Solving  these  equations  and  using  simple  geometry,  the  refraction 
compensation  transformat  ion  Is  defined  as  follows: 

r  ■  nnr  (22a) 

s 


H  •  $  — 


(22b) 


r 


where  (r  .  )  Is  the  transformed  polar  coordinate  representation  of  the 

s  s 

point  (r,  9).  i  is  the  look  angle  at  the  point  (r,  9),  Is  the  look 

angle  at  the  point  ( r  j ,  O) ,  and  i  is  some  arbitrary  constant.  The  value 

of  the  constant  i  depends  upon  the  application  of  the  transformat  Ion. 

For  example.  If  It  is  desired  that  the  transformed  rav  length  Is  equal 
to  the  true  rav  length,  then 


(1 


♦ 


/.dXrl 

/dx 


(23) 


where  the  line  Integrals  are  evaluated  along  the  rav.  Therefore,  the 
parameter  i  Is  dependent  upon  the  rav  direction,  or  look  angle,  $. 

To  analyze  the  validity  of  certain  approximat Ions  to  a  for  all  rays, 
two  atmospheric  regions  for  exponential  atmospheres  were  studied  - 
namely,  the  region  between  the  surface  and  1  kilometer,  which  can  be 
considered  as  fairlv  linear  In  nature,  and  the  region  between  the 
surface  and  10  kilometers,  where  the  exponential  nature  of  refractlvltv 
becomes  more  apparent. 


Recognizing  that 


j—  (nr)  -  n  (!♦!'), 

a  T 


(24) 


IS 


It  is  possible  to  approximate  »  for  rays  propagating  between  altitudes 
r |  and  r  ,  by 


a  -  f2  '  fl 
n2r2  "  nlrl 


Another  approximation  for  i  is  Riven  bv 
ri  ♦  r2  -1 

«  '  1 1  ♦  —  1  . 


where  (  and  ,  are  the  ray  factors  at  r^  and  r,,  respectively. 

Equation  (2b)  iias  errors  comparable  to  or  less  than  those  of  (2S)  for 
anftles  less  than  1".  but  (2S)  is  a  better  approx imat ion  for  angles  above 
1*,  especially  for  tin-  surface  to  10  kilometer  reRion.  A  comparison 
between  the  two  approx imat ions  is  Riven  in  Table  2.  It  should  be  noted 
that  the  errors  Riven  in  Table  2  are  for  unstratified  atmospheres,  such 
that  -0.  S  !'  0.0  (l.e.  ,  normal  propaR.it  ion) ,  where  the  refractivity 

can  be  approximated  either  bv  a  linear  or  bv  an  exponential  fit.  The 
larRer  errors  generally  occur  for  the  smaller  values  of  !;  this  implies 
that  extrers  care  should  be  taken  when  usinR  these  approximations  in 
superref ract InR  atmosphere*  unless  the  angle  is  significantly  large 
(i.c.,  larger  than  10*). 

Several  remarks  are  in  order  concerning  the  refraction  compensation 
t  rans format  Ion: 

a.  The  factor  of  n  in  the  radial  component  of  the  t ransf ormat ion 
is  crucial.  Although  n  Is  verv  close  to  unity,  the  gradient  of  n  la 

comparable  with  r  !  so  that  this  i ai  tor  is  Important  and  should  not  he 
de let  ed . 

b.  The  earth-center  angle,  •*,  does  not  appear  directly  In  the 
transformed  coordinate*.  Therefore,  extreme  care  should  be  taken  In  the 
use  of  the  transformed  earth-center  angle,  0^ .  This  difference  between 

and  Is  not  considered  in  the  effective  radius  or  effective  earth's 

M 

radiu'  transformations. 


. .  Tin-  radius  of  the  earth  does  not  enter  directlv  into  the 

calculations  if  the  transformed  coordinates.  The  parameter  r  equals 

r  ♦h  and  for  a  spherical  lv  stratified  atmosphere  is  the  radius  of 
e 

curvature  of  the  surface  of  constant  refractivity.  This  is  the  most 
serious  defect  in  the  effective  earth’s  radius  model.  The  same  relative 
errors  in  altitude  exist  between  the  effective  earth's  radius  and  the 
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refraction  compensat Ion  t  ransf ormat  Ions  ah  between  the  el  feet lve  earth's 
radius  and  tin*  effective  radius  t ransf orm.it  Ions  (l.e.,  approximately 


d.  Although  the  evaluation  of  the  parameter  i  does  depend  slightly 
upon  angle  if,  the  dependence  is  nowhere  near  as  great  an  Indicated  by 
the  effective  radius  t ran* format  Ion,  which  Is  Inaccurate  except  for  low 
angles. 


CONCLUSIONS 


The  effective  earth's  radius  model,  although  In  use  for  many  years. 
Is  essentially  In  error.  The  techniques  used  to  derive  It  do  not  produce 
a  t ransform.it  ion  that  maps  ravs  into  straight  lines.  As  a  result,  it  is 
possible  to  obtain  significant  errors  under  certain  applications.  The 
refraction  dispensation  transformation  does  produce  transformed  rays 
which  have  zero  curvature.  Also,  the  effects  of  the  transformation  on 
the  earth-center  angle,  which  are  Ignored  in  the  effective  earth's 
radios  model,  can  he  significant.  In  general,  the  author  feels  that  the 
use  of  an  effective  earth's  radius  contains  enough  errors  to  justify  Its 
removal  from  use  In  the  literature.  The  refraction  compensation  trans¬ 
formation  presented  in  this  paper  is  recommended  to  replace  the  effective 
earth's  radtus  t  rAnsform.it  ion  wherever  possible. 
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Appendix 

RADIUS  OF  THE  EARTH 


Generally,  the  earth  ran  b«  considered  as  an  oblate  spheroid  (i.e., 
a  sphere  flattened  at  the  poles).  Due  to  a  non-uniform  mass-distribution, 
the  direction  of  the  gravity  is  tilted  awav  from  the  center  of  the 
earth;  the  shape  of  the  earth  whose  surface  is  normal  to  the  direction 
of  gravity  is  called  a  gaoi d.  Over  the  land  the  radius  of  the  spheroid 
is  generally  less  than  the  radius  of  the  geoid,  while  over  the  sea  the 
opposite  is  generally  true.  Due  to  slight  irregularities  in  the  shape 
of  the  earth,  several  different  models  exist  for  the  sire  of  the  earth. 

The  four  values  usually  given  are: 

(a)  Equatorial  radius,  r  ; 

a 

(b)  Polar  radius,  r^; 

(c)  Mean  radius,  (2r  +  r.)/3, 

ii  n 

(d)  Elllpticlty  or  flattening,  e  •  (r^-r^)/^. 

The  more  commonly  occuring  values  for  the  above  values  are  given  in 
Table  AI .  The  various  models  are: 

(i)  '' irh  Spheroid  of  IRRS.  Computed  hv  the  English  geodesist, 

A.  R.  Clarke,  this  model  is  used  for  charts  of  North  America.  Since 
Clarke  did  not  clearlv  define  his  units,  the  D.S.  Coast  and  Geodetic 
Survey  adjusted  Clarke's  values  in  1880  by  adding  approximately  85  feet 
to  obtain  the  standard  values  given  in  Table  AI ,  which  are  used  for  all 
charts  and  maps  of  North  America. 

(il)  Clarke  Spheroid  of  18f>0.  These  values  are  new  estimates 
made  by  Clarke  in  1880;  however,  they  have  not  been  adopted  bv  the 
United  States. 

(ill)  International  Spheroid .  In  1909-10,  Hayford  (8)  con¬ 
ducted  measurement s  such  that  the  error  in  ie  equatorial  radius  is  ♦  18 

meters  and  the  error  In  c  1  is  ♦  0.5.  These  values  are  considered  to  be 
a  slight  Improvement  over  either  of  Clarke's  models.  However,  Helmet  (9] 
increased  these  errors  to  ♦  53  meters  and  ♦  1.2,  respectively. 

(iv)  5ir  Harold  Jeffreys  (1 952  and  1970).  Jeffreys  [101  is  a 
standard  reference  for  parameters  concerning  the  earth;  however,  the 
va 1  ties  given  therein  are  not  referenced  to  anv  experiment. 
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For  most  applications,  the  mean  radius  of  the  International 
Spheroid  Is  usually  adequate.  If  a  more  accurate  value  Is  necessary, 
for  a  latitude  (p,  assuming  an  elliptical  earth,  use 


r 

e 


2 

cos 


.  ,  2.2 

♦  r.  sin 
o 


* 


(A- 1 
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